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The relation between the angular diameter distance and redshift {dA-z relation) in a spherically 
symmetric dust-shell universe is studied. We have discovered that the relation agrees with that of 
an appropriate Friedmann-Lemaitre (FL) model if we set a "homogeneous" expansion law and a 
"homogeneous" averaged density field. This will support the averaging hypothesis that a universe 
looks like a FL model in spite of small-scale fluctuations of density field, if its averaged density field 
is homogeneous on large scales. We also study the connection of the proper mass of a shell with 
the mass of gravitationally bound objects. Combining this with the results of the dx-z relation, 
we discuss an impact of the local inhomogeneities on determination of the cosmological parameters 
through the observation of the locally inhomogeneous universe. 



I. INTRODUCTION 

The standard big bang model is based on the assumption of the homogeneous and isotropic distribution of matter 
and radiation. This assumption then leads to the Robertson- Walker (RW) spacetime geometry and the Friedmann- 
Lemaitre (FL) universe mode|| through the Einstein equation. This standard model has succeeded in explaining 
various important observational facts: Hubble's expansion law, the content of light elements and the existence of the 
cosmic microwave background radiation (CMBR) [||. 

The CMBR conversely gives a strong observational basis for the assumption of homogeneity and isotropy of our 
universe by its highly isotropic distribution together with the assumption that we are not in any special position in the 
universe (the Copernican principle). In fact, the deviation of our universe from a homogeneous and isotropic space is 
as small as ^ 10~^ at the stage of decoupling Thus our universe is well approximated by a FL model before this 
stage. On the other hand, the present universe is highly inhomogeneous on small scales; the density contrast against 
the cosmic background density is of the order of 10'^° for the sun, lO'^ on galactic scales, and of the order of unity even 
on the scale of superclusters. We have to go beyond FL models and linear perturbations in considering such systems. 

We usually regard that a FL model is a large-scale "average" of a locally inhomogeneous universe (averaging 
hypothesis). Even though the observational data are consistent with the picture that our universe is described 
well by a RW metric with small perturbations, we are not sure how to derive the background FL model from the 
inhomogeneous universe by any averaging procedure, or how the non-linear inhomogeneities on small scales affect 
large-scale behavior of the universe Q. Although one can derive a background FL model from observations of the 
nearby galaxies with any rule of averaging one likes, it is uncertain whether this background FL model (or its time 
evolution backward) agrees well with the highly homogeneous universe at the early regime. The discrepancies might 
appear, for example, in the estimate of the density of baryonic matter, the density parameter, the age of our universe, 
and so on. These still remain a non-trivial question to which we have to give a clear answer. 

Before proceeding to discuss these problems, wc should make the meaning of the term "average" definite. Averaging 
can be defined as a mapping between an inhomogeneous spacetime and a homogeneous one. This mapping is not 
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restricted to simple volume averaging; other averaging methods such as deformation of three-geometry according to 
the Ricci-Hamilton flow may be possible jij. By averaging, anyway, we expect that the large-scale (or coarse grained) 
behavior of the inhomogeneous spacetime is extracted. 

Averaging problem has been often studied from the viewpoint of the so-called back-reaction problem, i.e., how 
the small-scale inhoniogeneities affect the global dynamics when averaged on larger scales. Using the perturbation 
formalism and volume averaging, the back- reaction problem has been studied by several authors 1^-^. Apart from 
the problem whether the volume averaging is appropriate or not, one reason which makes the discussion of averaging 
unclear is that there exists no natural choice of time slice in an inhomogeneous universe. This leads to an ambiguity 
in the definition of averaged expansion rate or averaged density. Actually, there are even apparent discrepancies in 
the statements of the papers above. 

One possible prescription to avoid this difficulty is to treat observable quantities which we can give a clear definition. 
From this viewpoint, observational effects on the CMBR of small-scale inhomogeneities have been studied Jll| , p^ . 
These are significant in the sense that the back reactions on the observable quantities were discussed (see, e.g., llj] 
for a recent discussion of observable quantities in the Lemaitre-Tolman-Bondi spacetime). However, the problem how 
to determine the cosmological parameters by observing an inhomogeneous universe seems to have been overlooked. 

In Sugiura et aZ.(1998; hereafter Paper I) [Q, we investigated the distance- redshift relation in a spherically sym- 
metric dust-shell model, and compared it with that of a FL model. We focused on a highly (locally) inhomogeneous 
model; we prepared pressureless fluid distributed in discrete shells, which cannot be described by small perturbations 
of a FL background. We discussed the relation between the behavior of the dj^-z relation and the conditions on the 
initial time slice concerning the mass density and the expansion rate. We found that the distance-redshift relation 
observed at the center obeys a FL-like relation, even when there exist only several shells in the initial horizon scale, 
if the following conditions are satisfied on the initial spacelike hypersurface: the expansion law is homogeneous and 
the density which is averaged on larger scales than the inhomogeneity scale is scale- independent (i.e., large-scale ho- 
mogeneity of density field) . Here the averaged density is defined simply by dividing the mass of shells by the volume 
of the hypersurface. 

Our previous analysis was, however, limited to a spatially flat case (the Einstein de-Sitter model). There is a claim 
that the Einstein-de Sitter solution is in a special position as a solution of the Einstein equation; it is a kind of 
fixed point under renormalization group flow defined by the scaling properties of the Einstein equation |p^ . Thus, 
there is a possibility that our results were due to the special behavior of the perturbations of the Einstein-de Sitter 
model under an averaging procedure. In this paper, we treat non-flat cases and show that statements similar to 
those of our previous paper do hold. We also answer some questions unsolved in Paper I. In particular, the previous 
study of spatially flat cases could not answer which mass should be used in defining the average density, proper mass 
which appears in the stress-energy tensor of dust shell, or gravitational mass which specifies the parameter of the 
Schwarzschild spacetime; they are the same in the spatially flat model. Also it will be shown that the significance of 
the effect of inhomogeneities has curvature dependence. 

We here note three interesting points of the dust-shell model. First, its dynamics is exactly solved; it is not 
necessary to assume the existence of homogeneous background in order to obtain the evolution of matter distribution. 
Secondly, it treats a discrete mass distribution where the linear perturbation theory is invalid, and also can treat a 
highly general-relativistic situation where the scale of inhomogeneities are comparable to the horizon scale. In a sense, 
sphericall dust-shell model is a limit case of the Tolman-Bondi solution. We expect it can represent the Tolman-Bondi 
and the FL solutions when we take appropriate limits where the number density of dust shell goes to infinity, though 
no rigorous proof has been yet obtained. Thirdly, the motion of each shell can be extended even after shell-crossing 
occurs. This theme is discussed elsewhere |l5|,M. 




In order to get an insight into the inhomogeneous universe from the study of the dust-shell universe, we will also 
investigate the origin of the dust shell. Assuming that the dust shell is composed of highly bound objects, we consider 
the effect of the binding energy of each object on the dynamics of the universe. 

The organization of this paper is as follows. In the next section, we summarize the basic equations for the dynamics 
of a dust-shell universe, distance to the shells from the center, and redshift of the shells measured by an observer at 
the center. This is an extension of the treatment in Paper I. We give our results and discussion on d^-z relation and 
averaged density in Sec. HI. In this section, we also discuss the universe filled with gravitationally bound objects and 
the effect of the binding energy of each object on the dynamics of the universe. Finally Sec. IV is devoted to the 
summary. 

We follow the sign convention of the Riemann tensor and the metric tensor in |0| and adopt the unit of G = c = 1 . 




II. FORMULATION OF DUST-SHELL UNIVERSE 
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A. Equation of motion of dust shell 



First, we derive the expansion law of the dust-shcU universe. We consider a number of sphericaUy symmetric sheUs 
with a common center at r = (see Fig. 1). We label each shell by 1, 2, • • -,1, • • • from inside. The region enclosed 
by the first shell is the Minkowski spacetime and is labeled as the first region. Similarly, the region enclosed by 
the {i — l)th shell and ith shell , which is a vacuum, is labeled as the ith region. Each shell is infinitesimally thin, 
characterized by the surface stress-energy tensor S"''' which is given by 

S"'' = lim J T"'' dx, (2.1) 

where a; is a Gaussian coordinate (x = on the shell) in the direction normal to the trajectory of the shell. 

Since each region between the shells is a vacuum, the spacetime is described by the Schwarzschild geometry by 
Birkhoff's theorem. The line element in the ith region is written in the form 

ds^,=-(^l-^^dt^+(^l-^^ \r^+r^dn\ (2.2) 

where the parameter rrigi will be referred to as a gravitational mass (mgi — 0), and dO?' is the line element of a unit 
sphere. 

We derive the expansion law of a dust shell following p^-pl[|. Let rf" be a unit spacelike vector normal to the 
trajectory of the shell, and define the projection operator h^= 5\ — Ua'n!' ■ From the projected Einstein equation 

RabKh\ - 8^ [Tab - \9abT] h^^h^ (2.3) 



one obtains 



2- 



£nK,d + ^Rcd - KK,d = Svr ( T^bKh^ - ^h.^T ] , (2.4) 



where £„ is the Lie derivative along and '^Rcd is the three-dimensional Ricci tensor of the timelike hypersurface 
generated by the motion of the shell. The extrinsic curvature Kab is defined by Kab = —^h'^h'^£nhcd, and K ~ K^, 
T = T^. Integration of Eq. (E^) over an infinitesimal range along n° yields 



Kb - Kb = ( Sab - h^bS ] , (2.5) 



where the suffix denotes a quantity evaluated at the outside of the shell, and '— ' at the inside. Using Eq. ( |2.5| ) 
and the Gauss-Codazzi relation 2Gabn-°n'' = —^R + KabK"'^ — K^, one finds that the following relation holds for a 
dust shell: 

S-''{K++K-,)^0. (2.6) 

The stress-energy tensor of a dust shell is given by 

S"'' = (2.7) 

where ■u°_is t he 4 - velocity of the observer rest on the shell, and s is the surface energy density of the shell. Combining 
Eqs. (2_^), ( ^.6[) and (2/7), we obtain the following equation for the circumferential radius of the ith shell (the 
"expansion law" of the dust shell): 



dnV 2M, \ ( M^{i)V \ ml{i) 



dr J Ti I \ ms(i) y j Arf 
where M,; and M_ [i) are defined by 



(2.8) 



3 



M_(i) = mgi+i - mgj, 



(2.10) 



and r is the proper time which is measured by an observer rest on the sheU. The "proper mass" TOs(*) of the sheU is 
defined by 



(2.11) 



where Si = —S^{i). It can be shown this proper mass is a constant of motion by the conservation law, S'^.j, = 0, where 
the semicolon denotes the three-dimensional covariant derivative on the trajectory of the sheU. We shah assume the 
energy condition nis > 0. 

In this paper, we use a common proper time r for all the shells. The relation between r and the time coordinate, 
i, in the Schwarzschild spacetime is obtained as follows. First, note that two Schwarzschild time coordinates are 
assigned to each s hell: the it h sh ell has the time measured in the ith region and measured in the {i + l)th 
region. From Eqs.(2.5) and (2.6), we obtain 



dt 



(+)'■ 



dr 



— 2m, 



dt 



dT 



2m„ 



M_(i) _ m,(i) 
ms{i) 2ri 

M^{i) ^ m,(^) 
ms(i) 2ri 



(2.12) 
(2.13) 



The procedure to determine the origin of each time coordinate will be described later. 



B. Cosmological parameters and initial condition 



In order to specify a du st-shell universe, we have to fix the parameters in Eq. ( |2.8| ) and the initial hypersurface. 
We first rewrite Eq. ( |2.S| ) in the form corresponding to the Hubble equation of a FL model. We denote the initial 
circumferential radius of the ith shell by Xi, i.e.. 



on initial hypersurface. We define pi by 



(2.14) 
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M, = -■Kp^xC^, (2.15) 



and ki by 

M_(i) 
ms(i) 



E, = (2.16) 



Here Ei is the specific energy of the ith shell. Ei is positive through the space in the open and flat FL models and 
inside the maximum radius in the closed model. In this paper, we consider the cases where Ei > 0. 
Using these parameters, the expansion law of the dust shell can be written as 

We see that the first term behaves like a non-relativistic matter term in the Hubble equation of a FL model, the 
second and the third like a curvature and a radiation source term. From this point of view, pi and ki play roles of 
the "energy density" and the "curvature" , respectively. The radiation-likc term might be regarded as the effect of the 
binding energy of the shell |2^]. Further it is worthwhile to note that this radiation- like term is consistent with the 
result of Futamase |^ about the effect of the small-scale inhomogeneities on the global cosmic expansion. Seeing this, 
one may expect that the inhomogeneities tend to make the Hubble parameter larger compared with a homogeneous 
universe which has the same "energy density" of non-relativistic matter. However, this radiation-like term does not 
necessarily imply the larger Hubble parameter. In order to see the effect of this term on the Hubble parameter, we 
need to investigate the distance-redshift relation by solving the null geodesic equations and compare the result in the 
inhomogeneous spacetime and that of the FL model. We will discuss this point later. 



A dust-shell universe is specified if we set the parameters contained in Eq. (2.17), i.e., pi, ki, Xi, and an initial 



hypersurface. When we increase the number of the shells to infinity with pi and ki being finite and independent 
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of i (wc will mention this limit as "large N limit" ) , the dust-shcU universe approaches a FL universe if we take an 
appropriate initial hypersurface, as we will see in Sec. III. Then the parameters pi and ki agree with the ordinary 
energy density and curvature in the Hubble equation. 

We need to derive the expression of ti^±)i in terms of for the later use. The relations between t^^±^)^ and are 
given by [we denote ms{i) as m, in the remainder of this section] 



dt. 



(±)i 



dn {ri-2M,TE^mi),/Xi{ri) 



where 



Xiin) = {Ei^ -lyi"" +2Miri + 



4 ' 



2Mi 



and is the sign of dvi /dr. This is integrated to give 

U\± = m±{ri) + ii {2M, ± E^m^) \u 
where %± are integration constants. The functions Tj± (r) and G,± (r) are given for E'^ = 1 as 



2M, T E,mi) 



+ %±, 



Ti±{r) = 



3Mf 



MiT + 6M; ± -Mirrii - 



for £2 > 1 as 

Ti±{r) 



E,, 



Gi± (r) = { v'X^ + 2M, ± ^ }' , 

' E,{2Er - 3)i\/, ± {2Er - 1){E,- - 1)^ 



Ei^-1 



{Ef - 1)VW^ 



2Mi ± EiiTii 



Gi±ir) = 
and for £^ < 1 as 



h + {Ei^ - l)r + ^JE^1,/Xdr) 



1 TOi 



Ei 



tLi — 1 



■ E^i2E^-m,±i2E^^DiE^-l)^ i ^^^^^^ 

{E^-l),/l^ j 



r + k 



Gi±{r) = l^^{h-2MiTEimi){k + r) + ^{k + 2Mi±Eimi){h-r)^ , 

where h and —k are the roots of Xi{r) = 0; 

Mi + ^JMI + (1 - Ef)ml/A 



h = 
k 



Mi - ^Mf + (1 - Ef)m^/4 



(2.18) 
(2.19) 

(2.20) 

(2.21) 
(2.22) 

(2.23) 



(2.24) 

(2.25) 
(2.26) 

(2.27) 
(2.28) 



There is a coordinate singularity on the Killing horizon in the Schwarzschild coordinate; i(±)i becomes infinite on 
Vi = 2mgj. In order to avoid this coordinate singularity, and for further calculation, the Kruskal null coordinate is 
convenient since we will treat the null geodesies in this spacetime. Outside the horizon in the ith region, r > 2mgi, 
the Kruskal null coordinates are given by 



U = — 4w mgi(r — 2mgi) exp 



r - 1 



V = +4^mgi(r-2mgi) exp (^^^^ > 



(2.29) 
(2.30) 
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where U and V correspond to the retarded time and the advanced time, respectively. Using these Kruskal coordinates, 
the hne element in the ith region is expressed as 



dsj = -- 



cxp 



2to„ 



dUdV + r^dn^ 



(2.31) 



Similarly to the Schwarz schild time coordinate, two pairs of Kruskal null coordinates, ?7(±)i and V(±)i, are assigned 
to each shell. Using Eq. (2.2C), we obtain the Kruskal null coordinates labeling the ith shell in the form 



exp 



Ti +Ti±(ri) +%± 
AMi ± 2E,m, 



(2.32) 
(2.33) 



These coordinates, U(±)i and V(±)i, are finite on — 2mgi and are well defined also for < 2mgi. Hence, we will 
use the expressions (2.32) and (2.33) for any r^. For = +1, both J7(±)i and V(±)i are negative when 2mgi is larger 



than ri. In this case, the ith shell with < 2mgi is located in the white hole part of the Schwarzschild spacetime. 
On the contrary, for = —1, the ith shell with < 2mgi is located in the black hole part. 

The determination of the integration constants Ti± corresponds to the choice of the initial hypersurface. The 
procedure to construct the initial hypersurface we adopt is summarized as follows. First, we choose a unit spacelike 
vector which is directed outward (i.e., pointing the direction in which label i increases). Taking this vector as 
a starting tangential vector, we extend a spacelike geodesic curve until it reaches the second shell. This spacelike 
geodesic curve defines the simultaneous hypersurface in the region between the first shell and the second shell. Next 
we extend from this intersection towards the third shell another spacelike geodesic which starts with another spacelike 
vector at the second shell. This second spacelike geodesic generates a spacelike hypersurface in this region. Repeating 
this process, we complete the whole initial hypersurface. 

From the above procedure, the integration constant of Eq. ( 2.20| ) is determined as follows. In the ith region, the 
unit tangent vector, of the spacelike geodesic is written as 



1 



r = 1 1 



(2.34) 



and the other components vanish, where ej is an integration constant associated with the geodesic equation. Now a 
choice of initial time slice reduces to a choice of Ci. In Sec. Ill, we will see the connection between the choice of 
and the dA-z relation. Once we fix e^, the equation for the trajectory of the spacelike geodesic in the (t,r) plane is 
given by 



dt 



dr (r - 2mgj)v/(l + ef)r - 2m^ 
Integrating the above equation, we obtain 

where 



(2.35) 



(2.36) 



1 + ef 



+ ef)r - 2mgi + Ciy/r 



2mg,e,(3 + 2ej) ln(^(l + ef )r - 2mg, + ^(1 + ef )r) 
(l + e2)3/2 



(2.37) 



and Di is an integration constant. Initially, we set = t{-)i and i(±)i = 0. Using these relations, we obtain the 
integration constants, T{±)i, as 



i=2 



(xi - 2M, T E,m,) 



(2.38) 
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C. Redshift and diameter distance 



We consider a light ray which is emitted from each sheU toward an observer rest at the center. The light ray travels 
along a future directed ingoing radial null geodesic, where "ingoing" refers to the direction from a shell toward shells 
labeled by a smaller number. 

An ingoing radial null geodesic is specified by a constant coordinate value V. We denote the circumferential radius 
of the ith shell when it intersects the null geodesic as Ri. Labeling the outermost shell by M, the radius Rm is equal 
to a;M(thc initial radius of the Mth shell) and hence in the Afth region, V = V(_)jv/(2^m) is satisfied along the null 
geodesic. Thus, on the (M — l)th shell, the following relation holds: 

V(+)m-i{Rm-i) - V(^-)MixM)- (2.39) 

This equation determines Rm-i- By the same procedure, we obtain the circumferential radii of all the shells at the 
intersection with the null geodesic: 

= V(_),+i(i?,+i). (2.40) 

We can determine Ri from the given Ri+i through this equation. 

In order to derive the expression of redshift, we first write down the components of the null geodesic tangent in the 
ith region, k^{i), which is given in the Kruskal null coordinate as 

fc^(i) = — exp f-^) (2.41) 

and the other components vanish, where uoi is an integration constant associated with the geodesic equation. We 

require that the observed frequency of the photon at each shell is uniquely determined. The observed frequency, 
Wob(i)' the ith shell is given by 

c^ob(i) - -fc^(*)wf+)(«) = (2-42) 

Wob(« + 1) = -k^(iyi_^{i + 1) = ico,^^:^^^. (2.43) 
Eqs. (^l2|) and ( pl3|) gives the relation between ujoh{i) and Wob(* ~ 1), for i > 2, as 

0Joh{i) = f{i)iOoh{i-l): (2.44) 



where 



dVi_\i I dr 

For the first region, a direct calculation leads to 

c^ob(l) - + ^) ^ob(O) ^ /(l)^ob(O), (2.46) 

where a'ob(O) is the frequency of the light ray observed by an observer rest at the origin r = 0. Thus, using the above 
relations, we obtain the redshift of the light ray emitted from the ith shell toward the observer rest at r = in the 
form 

l + z(z)^^ = n/(,). (2.47) 

t^ob(O) ^\ 

Our next task is to find the angular diameter-distance c^a • The definition of dp^ is 

dA ^ J, (2.48) 

where D is the physical length of the source perpendicular to the line of sight, and 9 is the observed angular size. 
Since the space we are considering is spherically symmetric, the diameter distance from the observer at the center to 
the ith shell agrees with the circumferential radius Ri when the null geodesic intersects it; 

dAii)=R^. (2.49) 

Now we can calculate dA-z relation in the dust-shell universe using the relations ( |2lo| ), ( pl7| ) and ( ^l9| ). 
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III. RESULTS AND DISCUSSION 



A. Setting the parameters of dust-shell model 



As mentioned in Sec. Ill the choice of the parameters and the initial hypersurface determines the behavior of 
a dust-sheU universe. Since we are interested in cases which have a FL limit, we set the parameters so that they 
approach a FL model in the large N limit. 

We set the mass distribution of the shells as 



Pi = Pc (independent of i), 



and ki in Eq. (2.17) as 

ki — kc 

Using Pc and fcc, we define parameters i?shcii and rjj by 



shell 



„ Pc kc 



Further, we define 



Then, 



holl- 



kc^in- l)Hlcn- 



(3.1) 
(3.2) 

(3.3) 

(3.4) 
(3.5) 



In terms of FL models, -ffshcii, ru and Q may be regarded as the "Hubble constant", "Hubble horizon radius", and 
the "density parameter." 
For Xi, we put 



a;i(ri„it) = iAx, 



with a constant interval Aa 



where iVu is some positive integer. 



Aa; EE 



(3.6) 
(3.7) 



Before we proceed, we estimate the magnitude of the radiation-like term in Eq. (2.17). From Eqs. (2.£) and (2.15) 
with mgi = 0, we find 



n2(4n-3)r2rH, 



where n is a positive integer. From Eqs. (I2I0I) and (|2A|), we find ms{i) ~ (mgi+i — mgi)/Ei and thus obtain 



'ms{2n — 1) = 



ms{2n) 



Qn? — 6n + 1 



6n 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



Thus, when we consider a large N limit with fixing Xi — r-aii / Nn) , the proper mass Gms(i) = 0{N^'^) is regarded 
as a small quantity, compared with AI^{i) and nigi. That is, the radiation-like term of Eq. ( ^.17 ) is of order of 
the other terms, and can be neglected when A'h is sufficiently large. We note that when kc is positive (i.e., > 1), 
the factor fl/y^l — kixf in Eqs. ( ^.10 ) and ( [3.11 ) is larger than unity, and vice versa when kc is negative. This means 
that, in the dust-shell universe, the effect of the local inhomogeneities (caused by condensing masses to a shell) on 
the expansion rate appears larger in a closed model and smaller in an open model than that in a flat model which 
has the same iJsheii and Xi's. 
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B. Distance-redshift relation and averaged density 



As described in Sec.II, the choice of corresponds to choosing an initial hypersurface. We fix the expansion law 
to be homogeneous, try some choices of e^, and study which makes the distance-redshift relation behave like that 
of a FL model. We will consider which physical meaning is carried by that choice, especially in terms of averaged 
density. We have already analyzed spatially flat cases (i.e., kc — 0) in Paper I, and found that the distance-redshift 
relation agrees with that of a (spatially flat) FL model quite well when we choose so that the averaged density is 
homogeneous. The averaged density inside the ith shell p(i) is defined by dividing the "mass" contained within 
by the three- volume Vol(i) on the hypersurface up to r^. We may call the density field homogeneous when (0(i) is 
independent of i. 

We here write down explicitly the three-volume 



Vol(i) = 47r 



5 

T 2 



2m„ 



zdr 



207rm„,^ 



2m„ 



„5/2 



5^3/2^ 



3(1 + 6?) ' 6(l + e2)2 ■ 2(l + e2)3 



■In 



V(l + e?)r,-2mg, + V(l + ef)r 



(3.12) 
(3.13) 



For i = 1, Vol(l) is equal to Anrf/S. 

Let us study homogeneous density cases, as in Paper L Now we have two kinds of mass in defining the averaged 
density; proper mass ms(i) defined in Eq. (2.11) and gravitational mass rugi which determines the dynamics of its 
outer shells. 

We first take the hypersurface so that p{i) using the gravitational mass is homogeneous. We can define the 
gravitational mass of the jth shell to be the difference between the gravitational mass parameters of the neighboring 
regions: rugj+i — nigj [— Ejms{j)]. We add the gravitational mass of each shell up to i — 1, and add only half of the 
mass of the ith shell. The averaged density thus defined is given as 



J2 E,msU)\ /^Vol(j) = -{mgi + mgi+i) / ^ Vol(j 

j<i-l J j<i j<i 



(3.14) 



As one can see from the definition of pc, the choice of Vol(i) (and hence the choice of e^) which makes this averaged 
density homogeneous, i.e., independent of i (and equal to pc) is 



4-7r 

Vol(^) = y(rf 



(3.15) 



We plot the dA-z relation of the dust-shell models fixed in this way in Figs. 2-4 (choice A). The employed parameters 
are fl = 1.0, 0.9 and 1.1, iJsheii = 1-0, and the number of shells within the initial Hubble radius iVn is set to be 4 and 
10. The outermost shell M from which a light is emitted is chosen to be 2.5 x iVn- We also plot the dA-z relation in FL 
models, which is determined by the initial Hubble parameter Hi, initial density parameter fij, and the redshift of the 
initial time slice zi. These parameters are set using the parameters of the corresponding shell model as Hi — iJsheii, 
ill = ^, and zi equal to the redshift of the outermost shell for A^h = 10. In the flat case, the data of the dust-shell 
models agree with the FL relation as seen in Paper I. In the non-flat cases, however, the FL models chosen in this way 
do not approximate the dust-shell models. Even when we increase the number density of the shell, no improvement 
is obtained. Thus, using the gravitational mass of the shell in averaged density is inappropriate. 
Next, we try the choice using the proper mass in defining the averaged density, 



(3.16) 



After some manipulation, one finds that the volume element should satisfy the relation 



Vol(i) 



47r /rf 
3 



E^. 



Er 



(3.17) 
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in order to make the averaged density p{i) equal to pc (choice B). In a flat model, this choice is the same as choice 
A. The dA-z relations for open and closed models are plotted in Figs. 5 and 6. These plots show that the relations 
agree quite well with FL models. Thus, we should use the proper mass in de fining the a veraged density. We note that 



owing to the amplification of inhomogeneity which appeared in Eqs.( 3.1C ) and ( 3.11 ), closed universe shows slight 
deviation one can notice when compared with the other cases. 

For comparison, we try the orthogonal slice, since, in a FL model, the simultaneous hypersurface is orthogonal to 
the trajectory of matter (choice C). We choose Ci so that the vector is orthogonal to the trajectory of each shell. 
From the condition iau'^_^_){i — 1) = at r = r^-i, we obtain = dri-i/dr. Fig. 7 shows the dA-z relation in this 
choic^**]. We can see a mild deviation from the FL model. 

C. Discussion on the dust-shell model 

In paper I we concluded that the dA-z relation in a dust-shell universe behaves like a fiat FL universe, when the 
following conditions are satisfied: the expansion law resembles the flat FL model; the behavior of averaged density 
field is scale-independent when we increase the scale of averaging; the averaged density agrees with pc (defined by 
Eq. |3.1| ]). The discreteness nature of the shell model plays no significant role in our spherical model. 

However, our analysis was limited to spatially- fiat cases (Einstein de-Sitter model). In this paper, we have examined 
the non-flat cases, and confirmed that the above statement remains valid. Moreover, paper I could not answer which 
mass should be used in defining the average density, proper mass which appears in the energy-momentum of dust shell, 
or gravitational mass which specifies the parameter in Schwarzschild spacetime; the latter includes the gravitational 
potential energy and the kinetic energy of the shell. In the spatially fiat models, they balance and the two masses 
agree with each other. From the results of this paper, now it becomes clear that the proper mass should be used in 
defining the averaged density. This might at first sound strange since the geometry of each region is determined by 
the gravitational mass which includes the kinetic energy and the potential energy. Where has their information gone? 
Reexamining the Hubble equation of FL models and the expansion law of the dust-shell models, one notices that it 



is contained in the curvature term. In Eq. ( 2.16 ), the curvature term is expressed as the ratio of the proper mass to 
the difference of gravitational mass between the inside and the outside regions. Thus, the curvature gives the ratio of 
the total energy to the proper mass. 

In our model, the expansion law for the circumferential radius is homogeneous when there are enough number of 
shells; the effect on the expansion rate of density inhomogeneities is small (of order N^'^). Fixing the expansion law 
like FL, we studied the connection between the dA-z relation and the averaged density. We have found that if we 
make the averaged density (which is defined using the three- volume element of the hypersurface and the proper mass) 
homogeneous, the dA-z relation agrees well with that of a FL model. That is, there exists a strong connection among 
the homogeneous averaged density, the homogeneous expansion law, and the FL like distance-redshift relation, in 
spite of the discrete nature (local inhomogeneity of matter distribution) of the dust-shell model. This will support 
the "averaging hypothesis" that a universe behaves like a FL model in spite of small-scale fluctuations of density 
field, if its averaged density field is homogeneous on large scales. This implies that even if there exist large wall-like 
structures, our universe is approximated by a FL model on larger scales if the walls satisfy the above conditions. 

Now we give some comments on the small N cases; the cases where the radiation-like term in the expansion law 
is not negligible. One may expect that a FL model including radiation term can fit them, but this does not work. 
It should be noticed that the radiation-like term is inhomogeneous (i.e., dependent on i) in the cases considered in 
this paper. Then, there is no reason one may expect that the behavior of observables obeys a FL like relation; it is 
natural that we cannot fit them with FL models including radiation term. That is, it is impossible to approximate 
such an inhomogeneous model that has a significant large-scale inhomogeneity by a FL model. 

D. Implication on the universe filled with bound objects 

We can obtain further insight into the treatment of inhomogeneous universes by investigating the dust-shell universe 
from a different viewpoint, especially in connection with gravitationally bound objects. We discuss its implications 
on the estimation of the cosmological parameters of a locally inhomogeneous universe. 



** We only displayed an open model. There is no solution when > 1 which satisfies the conditions (3.6) and i^aW^+j (i — 1) = 
simultaneously for A'^h = 10, M = 25. 
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We start by considering the construction of a dust shell from small particles. In Appendix, we consider momentarily 
static initial data in which the intrinsic metric of the spacelike hypersurface is conformally flat and there are arbitrary 
number of compact objects. By construction of the solution, locations of the objects are arbitrary. Hence, by arranging 
an infinite number of infinitesimal objects with an appropriate procedure, we can construct an infinitesimally thin 
shell which is momentarily static. As discussed in Appendix, a spherical shell constructed by this procedure is likely 
to be regarded as a dust shell treated in this paper. Hence, by arranging sufficient number of such shells by the 
manner shown in this paper, we can obtain a system of compact objects which well imitates a FL universe. 

Next, let us consider a closed FL universe filled with baryonic dust fluid and a dust-shell universe (or a universe 
filled with compact objects) which have the same baryonic mass. The analysis of the dust-shell universe implies that 
the sum of the proper mass of all the shells in the closed dust-shell universe is the same as the baryonic mass of the 
closed FL universe, when the distance-redshift relations of these universes well agree with each other. On the other 
hand, by the investigation in Appendix, we find that the proper mass of a shell is the sum of the gravitational mass of 
the objects composing the shell. Here we should recall that the gravitational mass of an object is in general different 
from its baryonic mass; there is a gravitational mass defect and the difference between the gravitational mass and 
baryonic mass is recognized to be the binding energy of the object. Thus, if the shells are composed of highly bound 
objects, the sum of the gravitational mass of all the objects is much smaller than the baryonic mass. Therefore, even 
if the total baryonic mass of the universe composed of gravitationally bound objects is the same as that of the FL 
universe filled with the baryonic dust fluid, the distance-redshift relations of these universes can highly disagree with 
each other. 

Now let us study the difference between these universes with the same amount of baryonic mass quantitatively. 
The line element of the closed FL universe filled with the dust fluid is written as 

The solution of the Einstein equation can be written in the form 

a= — (1-COS77), (3.19) 

2M , , , , 

— {rt-sinr,), (3.20) 

where M is the total mass of the dust and < ^7 < 27r. In the dust-shell model, the mass M is given by the sum of 
the gravitational mass of all the objects in the universe if the shells are composed of bound objects; we denote it by 
Mq- On the other hand, in the FL universe filled with the baryonic dust fluid, M is the total baryonic mass Mb- 
When the two universes have the same amount of baryons, Mq < AI-q holds. 

Here we study how the Hubble parameter H and density parameter SI at a fixed age t = tq change when we change 
the mass parameter M. The Hubble parameter and the density parameter are defined by 

rr - 1 ^° ^ 37rsin?7 

~ adT 2M(1- cos 77)2' ' 

iira-^H'^ sm 77 

where M and 77 is connected through the relation tq = (77 — sin 77). Then the derivative of M with respect to 77 is 
given by 



dM 



37rTo(l-cos77) 

' 2(77-sin77)2 ■ ^^-^"^^ 



From Eqs.( |3.2lD , ( |3.22| ) and ( |3.23| ), we obtain 



dH 
'dM 

dn 

'dM 



Stt {3 (77 - sin 7;) - 77 (1 - cos 77) } 

2Af2(l_ cos 77)3 ' ^ > 

2(l-cos7;)(77-sin77) 
Msin^T] ■ 



One can confirm the positivity of dH/dM and hence the Hubble parameter increases with increasing mass M. On the 
other hand, the density parameter $7 is a decreasing function of M in the expanding phase while an increasing function 
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in the contracting phase. This imphes that if the dust shells are made of bound objects, the Hubble parameter of the 
dust-shell universe is smaller, and the density parameter is larger, than those of the FL universe filled with baryonic 
dust fluid (in the expanding phase) with the same amount of baryons and the same age, since the relation AIq < M-q 
holds. 

In the limit oi t/M — + (accordingly 77 ^ 0), the closed FL universe filled with dust fluid behaves like the 
Einstcin-de Sitter universe. The behavior of H^^dH/dM in this limit is easily obtained: 



1 dH 



HdM 



^ for ^^0. (3.26) 



30M M 



Using the above equation, the difference in the Hubble parameter of the dust-shell universe, i?DSj and that of the FL 
universe filled with the baryonic dust fluid, iJpL, is given by 

Thus, at the nearly Einstein-de Sitter stage, the effect of the mass defect of the objects in the universe is rather small. 
(The difference in the density parameters of these universes vanishes in this limit.) On the other hand, at the stage of 
the maximum expansion rj — n, the derivatives, H^^dH/dM and dil/dM, blow up. This is simply because H goes 
to at = TT and thus the diverging behavior itself has no serious consequence. However, this indicates the tendency 
that the effect of the binding energy of the compact objects becomes somewhat larger than the Einstein-de Sitter 
epoch when the curvature of the universe is comparable to the energy density of the dust fluid. 

The observation of CMBR strongly suggests that our universe was highly isotropic and homogeneous at least on 
the last scattering surface. Hence in the study of the universe in the early stages, the linear perturbation analysis 
is powerful. In order to perform the linear analysis, we need to fix the background FL universe, whose cosmological 
parameters are usually determined by the observation of our neighborhood. The universe observed today is, however, 
highly inhomogeneous and the inhomogeneities may prevent the correct determination of the background FL universe. 
The above discussion implies that taking account of the binding energy may be important in estimating the density 
parameter and the Hubble parameter near the maximum expansion, if the universe is filled with highly bound objects. 
On the other hand, when the universe is in the stage when the curvature of the universe is not dominant, the effect of the 
mass defect is rather small. We note that, however, in order to discuss the physical quantities of the inhomogeneous 
universe (e.g., the age of the universe) and its whole time evolution using a FL model constructed by the nearby 
observations, we have to know the behavior of the scale factor in the transition regime from the almost FL universe 
to the inhomogeneous one. This problem is now under investigation and will be given elsewhere [E3[. 



IV. SUMMARY 



We studied the behavior of d^-z relation in a spherically symmetric dust-shell universe where the mass distribution 
is discrete. Extending the treatment in our previous paper which only spatially flat models were considered, we 
analyzed non-flat cases. We compared the distance-redshift relation of dust-shell universe with that of FL models. In 
particular, we examined the behavior of the averaged density of the dust-shell universe when the two cIa-z relations 
are similar. We found that the d^-z relation observed at the center agrees quite well with that of a flat FL model if 
the following conditions are satisfied: (i) the expansion law of the circumferential radius of the shells resembles the 
Hubble equation of a spatially fiat FL model, (ii) the behavior of averaged density around the observer at the center 
is scale-independent as we increase the scale on which we take the average, and (iii) the averaged density agrees with 
the energy density of the FL model. In defining the averaged density, we take the total proper mass of the shells and 
divide it by the three-volume of the initial hypersurface. We noted that the choice of the initial hypersurface relates 
the expansion law to the averaged density. 

The effect of discreteness of mass distribution appears in the equation of motion of each dust shell. This effect 
becomes smaller as we increase the number density of shell, though we found that the positive curvature has tendency 
to enhance the inhomogeneity effect. We conclude that, in this spherical dust-shell model, the discrete nature of 
matter distribution plays no significant role in discussing the observed quantities such as c^a and z, as long as the 
expansion law and the averaged density field is sufficiently homogeneous in the sense described above. This supports 
the averaging hypothesis that a universe is described by a FL model if the universe is homogeneous when the density 
is averaged on a large scale than the scale of the inhomogeneities. This may also imply that even if there exist 
large wall-like structures, our universe is approximated by a FL model on larger scales if the walls satisfy the above 
conditions. 
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However, we have to keep in mind that our model is highly idealized and our analysis is limited only to spherically 
symmetric cases. In general, local inhomogeneities strongly affect the light propagation, giving rise to dispersions in 
the observed (Ia-z relation |24-26|. It will be also interesting to study cases when the light ray enters a shell in a 
non-radial direction. 

It still remains unclear whether three-dimensional discreteness have a significant effect on the dynamics of the 
universe. The study of the momentarily static initial data in Appendix strongly suggests that when the universe is 
filled with bound objects, the dynamics of the universe is determined by the gravitational mass density of the objects 
but not by the baryonic mass density. The gravitational mass of a compact object is in general different from its 
baryonic mass due to the gravitational mass defect. The effect of the mass defect on the dynamics of the universe is 
significant at the stage of the curvature-dominant phase while it is not significant at the early stage of the universe, 
i.e., the stage during which the universe behaves as the Einstein-de Sitter universe. However, in general, there are 
both the bound and unbound objects. Hence we should consider a situation including both objects and investigate 
their effects on the dynamics, which is left for our future work. 
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APPENDIX: MOMENTARILY STATIC INITIAL DATA OF A SPHERICAL SHELL COMPOSED OF 

BOUND OBJECTS 



In order to get an insight into the origin of the proper mass rus of a spherical shell, we consider momentarily 
static initial data of which the extrinsic curvature vanishes. The intrinsic metric of the three-dimensional spacelike 
hypersurface is assumed to be conformally flat, 

df = (Al) 

where a; is a position vector. Then the Hamiltonian constraint is written as 

A?A = -27rpHV'^ (A2) 

where A is the Laplacian operator in the flat space and pu is the energy density for an observer whose trajectory is 
normal to the spacelike hypersurface. 

We introduce the gravitational mass density deflned by 

PG = PhV''. (A3) 
Let us consider a situation where n "spherical" objects exist. We give pQ by 

n 

PG = Y.PGi^ (A4) 

7=1 

where 



PGI = PGi{\x-xj\)>0, (A5) 

for \x ~ xi\ < £i (radius of the Ith object), while pci vanishes for \x — xi\ > £i. The solution of the Hamiltonian 
constraint is then written in the form 

n 

ij = l + J2i^i, (A6) 

7=1 

where ^pI satisfies the equation 
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Ai^i = -27rpG/, (A7) 



with the boundary con dition, tpi for |.t — xi\ — > oo. 

Integrating Eq. (A7), the solution of the Hamiltonian constraint is easily obtained: 



i{x). (A8) 



POD n 

I=l''\S~x,\ Jo 

For the vacuum region, the above solution takes a simple form 
where the parameter m/ is defined by 

mi = An / dxx'^pGi{x). (AlO) 
Jo 

We also consider the proper mass m-pi of the /th object defined by 

mp/ = / d^xpnii^^ = / d^xpGii'- (All) 



If the compact objects are composed of dust fluid, mp/ is the conserved rest mass (baryonic mass). The above integral 
is easily performed to give 



,np, = m,(l + J, + l^^i^|, (A12) 



where 

Si = — / d^xpGi-ipi = / dzz^pGi{z) I dyy^'^ i dxx^pQi{x). (A13) 

mi J mi Jo Jo 

Now let us consider the "gravitational mass" of the /th object, which includes the gravitational binding energy. In 
order to obtain it, we replace the compact object by an Einstein-Rosen bridge with the same mass parameter m/. This 
means that we employ solution (A9) even inside the object. Now instead of the compact object, we have a "sheet" 
with an asymptotic region \x — xi\^ 0. In this asymptotic region, one can define the ADM mass (gravitational mass) 
Mi which corresponds to the total energy of the /th object. Brill and Lindquist showed that M/ is given by 

Mi^mi\l + ^y--^\. (AM) 



XI ~ Xj\ 



Then the gravitational binding energy E^^^^ of the /th object is naturally defined as 

Eii,,a = Mi-mpi = -mi6i. (A15) 

The above equation implies that 6i is the specific binding energy of the /th object. 

Here let us consider a situation in which compact objects with almost identical mass m/ ~ m are distributed 
homogeneously on spheres with a common center labeled hy i — 1, N. The ith sphere contains the objects labeled 
by / in the range rii + 1 < I < rii+i, where ni = 0, < n^+i, and njv+i is equal to the total number of the objects 
n. Note that the position vector of the /th object on the ith sphere satisfies \xi\ — Ri. Let us denote the sum of the 
parameter m/ inside the ith region as 

rii 

Mgi = ^mi. (A16) 

1=1 
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We consider a limit where the number of the objects on each shell goes to infinity with Mgi fixed: n^+i — — > cx) for 
all i's. By this procedure, we obtain a system of infinitely thin N spherical shells with a common center. This system 
has the same configuration as that treated in this paper (see Fig.l). In the limit where rii+i — Ui is very large, the 
mean separation, L, between the objects is given by 



while the mass parameter nij is 



L^J ^, (A17) 

y Ui+i - 



^ M^±l^. (A18) 

ni+i - n. 



We assume that the radius of each object is ij — ami, where a is a constant. If we consider an Einstein-Rosen bridge 
instead of the object, we choose the constant a to be 1/2. This assumption implies that in the large-(ni+i — rii) limit, 
ii/ L approaches to zero since it is proportional to (n^+i — ni)~'^/'^. Hence a shell obtained by this limit is extremely 
sparse. In order to know whether the shell obtained by this limiting procedure is a dust shell, we need to study the 
time evolution of this initial data. Although we do not investigate its time evolution here, it seems likely that a shell 
composed of infinite number of infinitesimally small objects is a dust shell by the following reason. Since the mean 
separation between the objects is infinitely larger than the radius of each object, a direct collision between the objects 
is impossible when the shell is shrinking. Further, since mj / L is extremely smaller than unity, the effect of the gravity 
of a nearby particle which causes a non-radial motion is likely to be very small. These two should work to keep the 
spherical symmetry and the radial motion of the shell during its evolution. 

In order to obtain the relation between Algi and the mass parameter rrigi of the Schwarzschild spacetime in the zth 



region, we investigate the line element in the ith region. The conformal factor in Eq. (Al) in the iih region is written 
as 

^ - A, + ^, (A19) 

where Ai is a constant and R= \x\. From the continuity of the conformal factor, we obtain a recurrent relation for 

v4.i as 



A, = v4,+i + (Mg,+i - Mg,) . (A20) 



Since the value of A^+i in the outermost region is unity, we can obtain Ai for 1 < i < A'" by the above equation. Eq. 
( |Aig| ) gives the line element as 



+ {dR' + R'dn') ■ (A21) 



2A,i?, 

Introducing a new radial coordinate X = AfR, the above line element becomes 

dp ^(^1 + ^ {dX^ + X^dn^) , (A22) 

which is the line element of the Schwarzschild spacetime in the isotropic coordinate. Then the mass parameter mgi is 
now trivially obtained as 

mgi = AiMgt. (A23) 

As reviewed in Sec. II, the dynamics of an infinitely thin she ll is treated by Israel's formalism. The equation for 
the circumferential radius ri of the zth shell is given by Eq. ( [2.8[ ): 

driV f M(i)V 2Mi m^(i) , 

The proper mass ms{i) of a dust shell is conserved during its time evolution. We will write down ms{i) usin g the 



mass parameters of the objects composing the shell, by comparing the above equation with the solution of (Al). The 
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momentarily static situation corresponds to the moment of maximum expansion dri/dr = 0. Hence from Eq. ( A24| ), 
we obtain the relation 



M_(i) 



2M, mlii) 



(A25) 



In the present situation, we know the gravitational mass rrigi and the circumferential radius r.; which is related with 
R^ by 



2i?. 



(A26) 



Hence Eq. (A25) is regarded as an algebraic equation to determine the proper mass ms(i) of the ith shell. The 
positive roots of this equation are given as 



ms(i) — m. 



1 



2me 



2mt 



(A27) 



where m^i < rngj+i is assumed. In order to make the meanings of the above roots clear, we consider Eq. ( ^.12 ). At 
the moment of the maximum expansion, ri should be larger than or equal to 2mgi+i. First we consider the case when 
ri > 2m^i+ i. T hen we easily find that dt(^^-^i/dT is negative for ms{i) = ms+ and positive for ms{i) — rris^. Together 
with Eq. (A26), the positivity of dt(^^^i/dT implies that ms(i) should be equal to mg- for Ri > mgi+i/2 (i.e., the 
direction of the time coordinate should agree with that of the proper time r in this asymptotic region), while it 
should be equal to mg- for Ri < mgi+i/2. When ri — 2ragi+i, i.e., Ri — mgi+i/2, the solution ms+ agrees with mg- 
and hence ms (i) i s equa l to m s± in this case. 
From Eqs.([A20|) and (|A23|), we find 



"^g. = + ^ (^g.+i - Mg,) \ Mg 



Using Eqs.( |A26| ), ( [A27| ) and the above equations, we obtain 



TOs(i) = (Mgi+1 - Mgi) Ai+ 



M, 



2i?. 



(A28) 
(A29) 

(A30) 



imposing dtf^^-^i/dr > 0. On the other hand, in the spherical-shell limit, the gravitational mass Mj in Eq. (A14) of 
the Ith object on the ith shell becomes 



Mi — > mi \^A,+i + 
Hence we find in the limit n^+i — — > oo for all the shells. 



i+l 



2R, 



rrisii). 



(A31) 



(A32) 



i=,ii+i 



This equation implies that the proper mass of the shell is the sum of the gravitational mass of each object, not the 
sum of the proper mass of each object. 
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(/+i)-th region 



i-th shell 

FIG. 1. Dust-shell universe. Between the shells the space is empty. 




FIG. 2. Angular diameter distance-redshift relation in spatially flat (fee = 0) dust-shell universe for choice A. Data points 
are shown by cross (+) and star(*). The number of shells within the initial Hubble horizon A^h is 4(-*r) and 10(+). The total 
number of shells M is taken to be 2.5 x A''h. The solid line shows the dx-z relation in a flat FL model with Hi = Hsheii- The 
redshift of the initial hypersurface is identified with the redshift of the outermost shell for the case A't = 25, i.e., zi = z{i = 25). 
As shown in Paper I, the deviation between the flat FL model and the dust-shell universe is almost unrecognizable. 
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FIG. 3. Same with Fig.l for open models (Q = 0.9). The sohd line shows the dA-z relation in an open FL model with 
n = 0.9, Hi = iJsheii, «i = z{i = 25). We see a large deviation. 
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FIG. 4. Same with Fig.l for closed models(f2 = 1.1). The solid line shows the dA-z relation in a closed FL model with 
Q, = 1.1, Hi = Hsheii, zi = z{i = 25). We see that this FL model does not approximate the dust-shell model. 



19 



12 3 4 

Z 

FIG. 5. Angular diameter distance-redshift relation in the open (SI — 0.9) dust-shell universe for choice B. Compare with 
Fig.3. The FL relation (solid line) agrees well with the relation of the dust-shell model. 
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FIG. 6. Angular diameter distance-redshift relation in the closed (f2 = 1.1) dust-shell universe for choice B. Compare with 
Fig.4. The FL relation (solid line) agrees well with the relation of the dust-shell model. 
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FIG. 7. Angular diameter distance-redshift relation in the dust-shell universe for choice C (open model with Q = 0.9). We 
see a mild deviation. 
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